Space-time phase transitions in the East model with a softened kinetic constraint 
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The East model has a dynamical phase transition between an active (fluid) and inactive (glass) 
state. We show that this phase transition generalizes to "softened" systems where constraint vi- 
olations are allowed with small but finite probabilities. Moreover, we show that the first order 
coexistence line separating the active and inactive phases terminates in a finite-temperature space- 
time critical point. Implications of these results for equilibrium dynamics are discussed. 



Glasses are formed when liquids are cooled below their 
melting temperatures. These supercooled liquids eventu- 
ally fall out of equilibrium and form glass [1, 2]. The 
transition between a flowing supercooled liquid and a 
rigid glass occurs with only a modest change in tem- 
perature. The resulting glass has a microscopic struc- 
ture which appears to be indistinguishable from its orig- 
inating supercooled liquid, despite large, but not abrupt, 
changes in its macroscopic properties. Several theories 
connect the glass transition to some kind of unusual equi- 
librium phase transition [3-6]. More recently, studies 
based on the dynamical facilitation theory of the glass 
transition [7, 8] have demonstrated connections between 
glassy behavior and dynamical phase transitions, which 
take place not in equilibrium ensembles but in space-time 
ensembles of trajectories [9-15]. Such phase transitions 
are not typically accompanied by equilibrium phase tran- 
sitions [9, 11], although some models do support both 
equilibrium and dynamical transitions [14, 15]. The dy- 
namical phase transitions are controlled by a biasing 
field, herein known as s, which couples to a measure of 
space-time activity, K . Here we investigate these phase 
transitions in a softened version of the East model [16]. 

Glasses arc often separated into two categories, labeled 
fragile and strong [17]. Strong glassformcrs arc those 
whose transport properties (such as relaxation time r or 
viscosity 77) follow an Arrhenius scaling with respect to 
inverse temperature, logr ~ 1/T. On the other hand, 
fragile glassformers are those that have transport prop- 
erties that depend on temperature in a more dramatic 
"super-Arrhenius" fashion. Within the dynamical fa- 
cilitation theory [7, 8], strong glass behavior is linked 
with diffusive relaxation processes while fragile behav- 
ior is usually correlated with hierarchical relaxation. In 
practice, glasses can exhibit a range of behaviors interpo- 
lating between strong and fragile behavior. In fact, it has 
been shown [18, 19] that relaxation-time data for glass- 
formers can be collapsed onto a single parabolic curve 
where logr ^ 1/T^. This collapse was theoretically pre- 
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dieted based on kinetically constrained models (KCMs) 
of glassformers with hierarchical relaxation [7] suggest- 
ing (at least within dynamical facilitiation theory) that 
most molecular glassformers exhibit some form of hier- 
archical relaxation. The softened East model has such 
hierarchical relaxation, which motivates its use in this 
study. In particular, we will compare the behavior of 
the softened East model with the softened Fredrickson- 
Andersen (FA) model considered in [20], which has dif- 
fusive relaxation mechanism and scaling reminiscent of 
strong glasses. Wc will discuss similarities between the 
dynamical phase transitions in these two models, as well 
as important differences that arise due to the combina- 
tion of softened kinetic constraints and hierarchical re- 
laxation with super-Arrhenius scaling. 



I. SOFTENED EAST MODEL 

The East model is a schematic lattice model for a frag- 
ile glassformer [16]. We work in one dimension but we 
expect our results to generalize trivially to higher dimen- 
sion. The thermodynamic properties of the model are 
those of a lattice gas [21] with N spins that take binary 
values rii = {0, 1} and an energy E = '^i' where 

J is the energy required to create an "excitation" (a site 
with TLi = 1). The spins do not represent the particles of 
the glass-former directly: instead they the describe the 
potential for mobility in a dynamically heterogenous ma- 
terial, after individual particles have been coarse-grained 
away [7]. "Excitations" are regions of space where re- 
laxation can (but may not necessarily) occur. Regions 
which locally cannot support relaxation are modeled by 
sites with rii = 0. 

An excitation on site i relaxes, i.e. rii = 1 — > = 
with rate [ri]o->.i = ACi, where Ci — e + rii-i is a facilita- 
tion function (or kinetic constraint [22]) that depends on 
the state of the left neighbor site, n^-i. The facilitation 
process is such that an excitation on site i — 1 facilitates 
motion on its "eastern" (right) neighbor, site i. On fixing 
the excitation energy J and the temperature T, the re- 
maining rates in the system are fixed by detailed balance: 
[rj]o->i = jXCi where 7 = e"-^/^. 

In what follows we set A = 1, which fixes our time unit. 
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and we take Boltzmann's constant fcs = 1 throughout. 
Having fixed our units in this way, the system at equi- 
Hbrium is described by two dimensionless parameters, 7 
and e. We imagine that violating the kinetic constraint 
requires an energy U, so that e = e~^/'^ . Hence, one 
may equivalently consider the system as a function of 
the reduced temperature T/ J and the ratio of energies 
U/J. 

In terms of the equihbrium dynamics on decreasing the 
temperature (at fixed U and J), the soft East model ex- 
hibits two crossovers. At high temperature, facilitation 
effects are irrelevant, as usual in KCMs [22], with facil- 
itation becoming relevant around the onset temperature 
To ~ J. However, in the presence of a soft constraint, 
facilitation effects are also irrelevant at very low temper- 
atures, as we now discuss. The facilitated dynamics of 
the "hard" East model (with e = 0) lead to a time scale 
'''fac(^) for dynamical relaxation on a length scale t. the 
relevant rate is 



II. ACTIVE AND INACTIVE SPACE-TIME 
PHASES 

In order to investigate dynamical (space-time) phase 
transitions in the soft East model, we define an or- 
der parameter K that measures the activity in the sys- 
tem [9, 11; 26-29]. The activity is the total num- 
ber of configuration changes in a trajectory, K = 

Y^iLi Y.t=o ("-i(^) ~ + where tobs is the length 

("observation time") of the trajectory, and St is an el- 
ementary time step. Since we work in continuous time 
[30] we take the limit where St — s- 0. 

The activity K distinguishes between active and inac- 
tive phases (the active phase is characterized by a large 
value of K while the inactive phase is characterized by 
a small value of K). While K is an order parameter 
in space time, it is analogous to an equilibrium order pa- 
rameter such as density, which distinguishes between two 



rfac(^)-'^exp(-(J/r)log2^), 



(1) 



which is valid as long as ^ < c'^^ and J/T <^ 1. The 
equilibrium relaxation time of the model with e = scales 
as TEast ~ e'^ / ^^"^ [23] , comparable with the predic- 
tion [24] obtained by setting £ equal to the typical spacing 
between excitations, e"'/-^. When e > 0, one expects an 
additional relaxation mechanism via constraint violation, 
with rate t^^^ e exp(— J7/T). Thus, the structural 
relaxation time r in the soft East model may be esti- 
mated via 
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At low temperatures, the facilitated time scale diverges 
faster than Tgoft, and the softened mechanism dominates. 
(This is opposite to the behavior in the soft-FA model, 
where the facilitated relaxation dominates as T — >■ 0, as 
long as J > 3U [20].) The crossover to the softened mech- 
anism cuts off the super- Arrhenius divergence of t in fa- 
vor of an Arrhenius law, reminiscent of strong glasses. 
This crossover takes place at 



J2 



2[/ln2 
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and is accompanied by a saturation of the equilibrium 
dynamical length scale at ^sat ~ 2'^/''. However, in con- 
trast to the crossovers from fragile to strong behavior in 
other KCMs [7, 25]), we note that the activation energy 
for relaxation here increases monotonically with decreas- 
ing temperature. In experiments on glass-formers, we 
expect ?7 to be a large energy scale, so that is likely 
to be beyond the limit of supercooling of the liquid. How- 
ever, this low temperature crossover will be relevant for 
our discussion of phase transitions in the soft-East model 
below. 
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FIG. 1: (A) Phase diagram for the Id soft East model in the 
(s, e) plane at J/T = 0.75 (so 7 = 0.47). The solid line is the 
phase boundary between the active and inactive phase. The 
dashed line is the continuation of the symmetry line (4) into 
the 1 phase region. The red circles indicate a state point with 
e = 5 ■ 10~* < ec, on the coexistence line, in the two phase 
region. The blue triangles indicates a simulation point where 
e = 5 • 10~^ > £c- The black x indicates a state near the 
critical point: e — 1.5- 10^'' ~ ec. (The precise location of the 
critical point is not known for this model). The inset to (A) 
shows histograms of the intensive activity k for the three state 
points in (A). (B) Plots of average intensive activity (fc(s)) as 
a function of field s, for the values of e given in (A). 
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FIG. 2: Sample trajectories from the three state points identified in Fig. 1, taken from near the centers of the distributions 
P{k) {k ~ 0.1). Thus, for t < tc these trajectories are rare, coming from the trough in the histogram that lies between the 
two stable basins. Active sites are colored (rii = 1) and inactive sites are white {jii = 0). (A) Trajectory with e < e^. showing 
space-time phase separation. (B) Trajectory at e ~ tc where the phases are still identifiable but the clusters no longer have a 
sharp interface. (C) Trajectory at e > ec showing a single homogeneous phase. 



phases such as liquid and gas. However, while the equi- 
hbrium ensemble for a liquid-to-gas phase transition is a 
set of configurations, the ensembles we consider here are 
collections of space-time trajectories. (Mathematically, 
we are considering the "large deviations" of the order 
parameter K [27, 31].) 

To access space-time phase transitions, we introduce 
an intensive biasing field s which couples to K. This 
defines a nonequilibrium ensemble of trajectories [9- 
11, 26, 27, 32] known as the s-ensemble [31]. For an ob- 
servable A, we use {A) g to denote the expectation value 
of A in the presence of a biasing field of strength s. Sim- 
ilarly, {A)q denotes the equilibrium expectation value. 
These expectation values are related through 
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Z{s,tohs) 



(3) 



where .Z(s,tobs) = (exp(— sA'))o is the partition function 
for the s ensemble. 

While the field s has no direct physical interpretation, 
the formalism we use implies an equivalence between en- 
sembles with fixed K and those with fixed s [31]. That 
is, the field s acts in the same way as a constraint on K. 



III. COMPUTATIONAL SAMPLING OF 
SPACE- TIME PHASES 

To harvest ensembles of trajectories we use transition 
path sampling (TPS) [33], as in [12, 20]. We employ two 
basic move types to generate trajectories from an initial 
trajectory: shooting and shifting in both the forwards 
and backwards directions. New trajectories are accepted 
with a probability proportional to exp(— siiT) - analogous 
to accepting and rejecting configurations with probabil- 
ity proportional to exp(— in standard configurational 
Monte Carlo dynamics. 



For the softened East model, the s-ensemble depends 
on only three parameters: 7, e, and s. Fig. 1(a) shows 
a space-time phase diagram in the (e, s) plane, for fixed 
7 (specifically, J/T 0.75, so 7 = 0.47). On the s = 
axis, the system undergoes equilibrium dynamics, while 
the e = axis is the 'hard' East model. The solid line 
on the phase digram indicates phase coexistence between 
inactive and active phases, as estimated from our numer- 
ical results. On increasing s, one crosses the phase co- 
existence line at s = s* and the activity in the system 
changes discontinuously from an "active" state to an in- 
active one: this is accompanied by a jump in the value 
of A' [see Fig. 1(b)]. 

To characterize these phase transitions in the s- 
ensemble, we perform finite-size scaling, by varying the 
system size and the observation time tobs. As in [20], 
there are two analytic results that improve the quality 
of our finite-size scaling analysis. The theoretical anal- 
ysis is given in Sec. V: here we simply state the results 
and explain how we use them in conjunction with our 
TPS scheme. Firstly, we perform our TPS simulations 
at s = s*, which allows accurate estimation of the prop- 
erties of the phase transition. The transition point s* 
can be located exactly because there is a hidden symme- 
try of the s-ensemble for the soft East model, similar to 
that found in the soft-FA case [20]. This implies that s* 
satisfies 



1 + 7 
2e + 1 



= v/(l - 7)^ + 476-2"* 



(4) 



A second result that aids finite-size scaling analysis in- 
volves the boundary conditions used in the s-ensemble: 
as in [20], it is natural to take periodic boundary condi- 
tions for the N spins in the system, but the boundary 
conditions in time are not periodic, with the initial and 
final configurations of the trajectory being free to fiuc- 
tuate independently. This free boundary favors activity 
near initial and final parts of the trajectory [31], which 
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frustrates convergence of the limit of large <obs ■ To coun- 
teract this effect, we bias the initial and final conditions 
of the trajectory, as follows. We define an angle a by 

tana = — (5) 

1-7 

with < a < 7r/2, and wc define 

5sEast = ln[tan(Q;*/4)V7] • (6) 

where a* is the value of a when s = s*. We then in- 
troduce an extra bias on the s-ensemble, which depends 
on the total number of excitations in the initial and final 
conditions of the trajectory: in the resulting ensemble 
then the average of observable A is 

(yle -^-^+9 =Ea=t [A^(0) +AA(t„b= )] ^ 
s,symm " 'V { 4- \ 

(7) 

^symV^; ^'obs ) 

where M(i) ~ X^i^i ^A^) the number of excitations at 
time i, and Zsy,„(s, Us) = (e-"-^'+9=E„t[AA(o)+A^(tob3)])g is 
a normalization factor. As discussed in [20], this extra 
bias ensures that the s-ensemble at s = s* is fully sym- 
metric between the active and inactive phases, and we 
term it the symmetrized s-ensemble. 

In the limit where the observation time iobs is infi- 
nite, this symmetrized ensemble is equivalent to the en- 
semble without boundary constraints, (^)s,symin ^ 
as tobs CO [20], as long as observable A is not domi- 
nated by the initial and final parts of the trajectory. This 
boundary biasing condition is straightforwardly included 
in the transition path sampling algorithm. 

IV. NUMERICAL RESULTS 

The order parameter for the phase transitions we con- 
sider is the intensive average activity (fc(s)), 

(/c(s)) s . (8) 

J V lobs 

As discussed in the previous Section, our main results are 
summarized in Fig. 1 where we show sample numerical re- 
sults for fc(s) (panel b), and the phase diagram (panel a) 
that we have obtained by finite-size scaling of the behav- 
ior of k{s). In particular, there is a first-order space-time 
phase transition represented by a solid line, which ends 
in a critical point as (sc, ec). That is, the phase transition 
that is known to be present as e = [11] still exists in 
the presence of soft constraints, as long as e < ec. How- 
ever, when the constraints are too soft (e > Cc), the phase 
transition disappears and the system shows a smooth re- 
sponse to the field s, and no phase transition. The exact 
value of tc is not known for the soft East model, but 
we have bracketed its location by identifying two state 
points, one where the phase transition occurs, and an- 
other where the response to s is smooth (non-singular). 



The inset to Fig. 1(a) shows probability distributions 
P(fc), where fc = K/Ntohs is the (intensive) activity per 
space- time volume. The histograms correspond to sim- 
ulations performed at the conditions highlighted by the 
symbols along the symmetry line. These were chosen to 
lie in the two phase region far from criticality, near crit- 
icality, and in the one phase region which does not sup- 
port a phase transition. In the first order region e < Ec, 
there arc two distinct peaks corresponding to the inac- 
tive and active phases. Near the critical point e w ec the 
two peaks are broadened and very flat. In the one phase 
region e > there is one distinct peak. 

Figure 1(b) shows the jump in {k{s)) as the field s is 
varied. It is useful to compare these results with those 
expected for an equilibrium phase transition in a ferro- 
magnetic system (see also [9, 10]). The behavior of fc(s) 
is analogous to the change in the magnetization, M, of 
a ferromagnet as the field strength, h, is changed. Thus, 
e < Ec corresponds to T < Tc in a ferromagnet and shows 
a jump in the order parameter (from a high k active phase 
to a low k inactive phase). On the other hand, e ~ Cc 
corresponds to T ~ Tc in a ferromagnet: it shows a steep 
crossover between the phases (for the system sizes con- 
sidered here). 

In Fig. 2, we show three trajectories which have k ~ 
0.1, approximately corresponding with the center of the 
distributions. These trajectories are harvested under the 
same conditions as the marked points in Fig. 1(a). In 
Fig. 2(a) a trajectory for e < Ec is shown: it may be seen 
from the inset to Fig 1(b) that this is a rare trajectory, 
coming from near the minimum of P{k). There is a clear 
segregation between the two phases as the space-time 
"surface tension" [20] is large. In Fig. 2(b) a near-critical 
trajectory with e « Ec is shown. Here, there is still dis- 
tinct phase segregation but the boundaries between the 
phases have become amorphous, corresponding to the di- 
minished "surface tension" near criticality. In Fig. 2(c) 
we show a trajectory with e > ec. Here, there is no phase 
separation since the system in a one-phase region of the 
phase diagram. 

As usual for phase transitions in space-time [31], 
the first order jumps in fc(s) become singular only as 
N,to\ys — >■ oo. For finite systems undergoing a first-order 
phase transition, the susceptibility x = —d{k{s))/ds eval- 
uated at s* is proportional to the system size, iVtobs [34]. 
This scaling is shown in Fig. 3. In Fig. 3(a) we show the 
finite size scaling for the first order transition for e < Ec- 
As the system size increases, the transition sharpens. In 
Fig. 3(b) we show the scaling near the critical point. 
Here, the crossover sharpens as iViobs increases, but there 
is no clear jump in (fc(s)). In the inset to Fig. 3(b) we 
plot x*j the value of the susceptibility evaluated at the 
transition point s* as function of system size Ntoi^s- For 
e < ec and e w ec the scaling grows linearly. For e > e* 
there is no increase with system size. 
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V. THEORETICAL ANALYSIS 
A. The master equation and its symmetries 

The soft East model may be studied analytically fol- 
lowing the methods used for the soft FA model in [20]. 
We now use these methods to derive the relations (4) and 
(6) that we used to facilitate our numerical studies. To 
do so, we start with the master equation 

dtP{C,t)^~r{C)P{C,t)+ W{C' ^C)P{C',t) 

C'(#C) 

(9) 

where P(C,t) is the probability of observing a config- 
uration C of the soft east system at time t. Also, 
^(C) = J2c' — > C) is the exit rate to leave the 
current configuration C, sums over C run over all pos- 
sible configurations of the model, and W{C' C) is the 
transition rate from C to C. 

This master equation can be represented compactly 
using a spin-half representation of the system's config- 



e < Cc 

tobs = 640, TV = 60 + 

tobs= 1280, 7V= 60 X 

tobs = 2560, TV = 60 * 

tobs = 51 20, AT = 60 □ 
Us = 640, TV =120 

Us = 1280, TV = 120 A 



0.001 0.002 0.003 0.004 

(b) e ^ ec 





© £ < 












A e > 








y 






••k 





LS-IO'^ 3-10=^ 

tohs -N 



0.004 0.008 0.012 

s 



FIG. 3: Finite size scaling of the intensive activity (fc(s)). 
In panel (a) we show e = 5 ■ 10~* < tc and in (b) e = 
1.5 • 10~^ ~ Ec. Inset to (B) shows scaling of the suscep- 
tibility X* = ~d{k) /ds\^, for various system sizes iobs • A*', 
with symbols and colors corresponding to thee values shown 
in Fig. 1. The dashed lines indicates linear scaling associated 
with a first order phase transition (for e < ec and the constant 
value of susceptibility expected in a single phase (e > £c). 



uration space. The ground state (where = for all 
i) is denoted by jfi), with the site variables {n-i} now 
represented by N spin-half variables. Hence, all possible 
configurations {n-i} can be represented via Pauli matri- 
ces CT^'^'^ acting on \Q) [35], the raising and lowering 
a~ operators arc defined as = \{<j'f ± crj') Thus, a 
configuration {ji;} can be represented as 



N 



|{nj)=^(a+r[0) 



(10) 



We then write \P{t)) ^ J2c which allows us to 

rewrite the master equation (9) as 



|[P(t))=W[P(t)) 



(11) 



where W is a linear operator whose matrix elements are 
the transition rates, W{C — s- C) and escape rates r(C). 
For the soft East model this operator is 



= ^(n,_i + e)[(l - a+)ar + 7(1 - ar)a^ 



(12) 



This master operator is similar to that of the soft-FA 
model given in Ref. [20]. However it does not contain 
the symmetrized portion ( i o j ) as the East model is 
inherently asymmetric, nor does it contain the diffusive 
term that was introduced to the soft-FA model: introduc- 
ing such a term in the soft-East model does not simplify 
the analysis in the way that it did for the soft-FA model 
so we do not consider it here. 

To analyse the soft East model in the s-ensemble, we 
follow [31]. First the probability distribution P{C,K,t) 
is defined as the probability of being in configuration C 
at time t, having already accumulated an activity K in 
the time between time and t. The probability to be 
in configuration C at time t at a field strength s is then 
P(C, s, t), which can then be written as a reweighted sum 
over over all possible accumulated K values: 



P{C,s,t) = Y,PiC,K,t)t 



-sK 



(13) 



K 



The equation of motion for |P(s, t)) ^ P{C, s, t)\C) is 
d 



dt 



\P{s,t)) =Wis)\P{s,t)) , 



(14) 



where now the matrix elements of W(s) include the bi- 
asing field s, through 

W(s) = 5](n,_i + e)[(e-^ - a+)ar + l{e-^ - ar)at] . 

(15) 

This resembles Eq. (12), with the additional factors of 
modifying the rate of changes of state. 
This representation of the model allows us to derive 
the symmetry condition (4), following [20, 36]. Writing 
the energy operator E ~ J^i'>ii (recall 7 = e~'-'^'^), 
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we symmetrize W(s) using H(s) = c^^'^'^W{s)e'^/^'^ , so 
that 

i 

A + 7 A 

X ( — h (Ji cos a — erf sin a j (16) 

where = + 7)^ — 47(1 — e^^'*) and a was de- 

fined in (5): we have (cos a, sin a) = 77(1 — 1 ^'^'^^ ^ ■ 
Next, making a (site- independent) rotation of the 
quantum spin matrices, (cr^ , cr*' , cr^ ) — > (—cr^ cos a — 
(7^ sina, — (T^, CT^ cosa — cr^sina), one finds H(s) ^■ 
EI-w(s) with 

IHIw(s) = X/ (-"^ '''^'-i cosa - crfLi sina) 

x(i±2.<,r) (17) 

Comparison with (16) shows that (17) is a symmetrized 
master operator for a "West model" , similar to the East 
model except that spin i facilitates motion on its left 
(west) neighbor, spin i— 1. This West model has the same 
parameters as the original East model if 1 -I- 2e = ^i^: 
this corresponds to the symmetry condition (4). The spin 
rotation and the symmetrization of W(s) preserve the 
eigenspectrum of this operator, meaning that the model 
is 'self-dual', with a symmetry that may be spontaneously 
broken. Formally, there is a similarity transformation 

Hw(s) = f/lH(s)[/„ (18) 

where = Uji-e'""^ ^'^'^jY- see also [20, 36]. When the 
West model described by IHIw(s) has the same parameters 
as the East model described by IHI(s), this represents a 
symmetry of the model, which is broken spontaneously 
at the spacetime phase transitions. 

To understand the boundary bias introduced in (7), 
it is useful to write Z{s,tohs) ~ (e|e'**(^^*°''=|eq), with 

(e| - (0| n,(l + Icq) = n, We also 

introduce the general coherent state \(f>) = J|^(cos(/) -t- 

(T^ sin0)|ri), parameterized by the angle (f>. Then, sym- 
metrizing W, one arrives at 

Z(s,tob.) = (</>o|e^(^)*°^=|0o> (19) 
with tan^o = ^/j and < 00 < tt/2. From (18), one has 

Z(s,tobs) - ((a/2) - 0o|e^*(^)*-=|(«/2) - M (20) 

where we used U^\(j)) — | (a/2) — 0), which follows from 
the definitions of Ua and 

The key point is that the partition function (19) is fully 
symmetric under this similarity transformation only if 
two conditions are met: (i) ]HIw(s) should describe a West 
model with the same parameters as H(s) and (ii) the 



boundary terms are invariant under the transformation: 
(a/2) — 00 = </'0i so that (20) coincides with (19), up to 
the replacement of E[w(s) by H(s). It is easily verified 
that while (4) ensures condition (i), it does not simulta- 
neously ensure condition (ii). To simultaneously satisfy 
both conditions, we write 

^sym(s,tobs) (e 1 6^ e*^(')*°''= 6^'='=-' |eq) 

(21) 

Symmetrizing W(s) and using the expression (6) for 
ffsEast, one finds 

^sym(s,tobs) = y-^(a/4|e«(^)*°-|a/4) (22) 

with y = (1 + 7)cos^(a/4). Applying (18) yields 
-^sym(s,^obs) = c"^ (a/4| e"'w('')*ob= |q,/4^ ^ indicating that 
the partition function is fully symmetric under the sim- 
ilarity transformation, as long as (4) holds. Physically, 
this means that the initial and final conditions in the 
symmetrized s-ensemble have no additional preference 
for the active phase over the inactive one, ensuring that 
the probabilities of active and inactive states are equal 
within the ensemble that we sample. This facilitates ac- 
curate numerical characterization of the phase transitions 
in the system. 

B. Scaling of s* and ec with temperature 

As well as exact symmetries of the master operator, 
we also identify low-temperature scaling properties of the 
softened East model. For low temperatures, 7^1, the 
symmetry condition (4) reduces to s* w 2e/7. As pro- 
posed in [37], this scaling of s* with system parameters 
may be obtained by comparing the difference Ak in ac- 
tivity between the two phases, and the escape rate 70 
from the inactive phase. Specifically, we expect 

s* « 7o/Afc. (23) 

At low temperatures, the activity difference between the 
phases is given simply by the equilibrium activity per site, 
so Ak = 27^ A. The escape rate (per site) is 70 ~ eX^, 
since a single unfacilitated spin flip in a large inactive 
region typically leads to relaxation of the entire inactive 
region. Thus, the prediction (23) works well in the soft- 
ened East model. (The same analysis also applies in the 
softened FA model, except that Ak is larger by a factor of 
2, so the coexistence field s* is commensurately smaller.) 

It is instructive to generalize the mean-field analysis of 
the soft-FA model [20] to the East model. In both soft-FA 
and soft-East models, the mean- field analysis gives the 
correct low-temperature scaling for s* . However, at this 
mean field level, one finds that the critical point scales as 
Sc ~ 0(1) and Ec 7, as 7 — > 0. This analysis seems too 
simplistic, especially given that the mean-field analysis 
does not capture the hierarchical relaxation in the East 
model, and the mean-field predictions do not seem con- 
sistent with the small values of Ec and Sc shown in Fig. 1. 
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FIG. 4: Schematic phase diagrams in the (s, T) plane, (a) The 
"hard" East model {U — >■ co, varying T/J and s). The system 
is in the inactive phase for all s > and in the active phase 
for s < 0. The equilibrium dynamics are characterized by 
diverging length and time scales as T — 0. (b) The softened 
East model with fixed U / J > 1, varying T/J and s. For 
an intermediate range of temperature, there is a first order 
phase transition at s > 0. This is also the temperature range 
in which facilitation effects are strong. Facilitation effects 
are weak both for large T (above the onset temperature) and 
for small T (below Tx): there is no phase transition in these 
weak-facilitation regimes. 



A more appropriate estimate of Cc can be obtained from 
dynamical scaling at equilibrium in the East model: scal- 
ing requires that requires that Tgoft/TEast should remain 
constant as 7 is reduced. Assuming that the large de- 
viations of K have the same scaling as the equilibrium 
fluctuations, one therefore expects 



exp 



(log 7)^ 
21n2 



(24) 



with Sc as constrained by the symmetry condi- 

tion (4). These predictions indicate that the existence of 
a space-time phase transition requires that the probabil- 
ity of constraint violation must remain extremely small 
at low temperature, at least in this model. 



VI. IMPLICATIONS OF THESE RESULTS FOR 
THE GLASS TRANSITION 

While there is no direct experimental realization of the 
s-field, we believe the results presented here have impli- 
cations to real glassformers. Again, we are motivated via 
analogy to the ferromagnetic system, as discussed in [10]. 
Imagine a ferromagnet [21] at temperature T < Tc and 
in a magnetic field h > 0. The dominant phase is or- 
dered and aligns with the field h. However, there remain 
fluctuations which arc aligned opposite to the preferred 
phase. The size of these fluctuations and their shape are 
controlled by the free energy difference and surface ten- 
sion between the two phases [21]. Near Tc, the character 
of these fluctuations is influenced by the nearby critical 
point and the fluctuations become larger and their shape 
becomes less well defined. 

For 'hard' KCMs (with e = 0), the coexistence be- 
tween the active and inactive phases occurs exactly at 
s* =0 [11]. On the other hand, for the softened models 
presented here and in Ref. [20], the coexistence moves 



towards s* > 0. The equilibrium dynamics of the sys- 
tem corresponds to the active phase which dominates the 
system. However, for equilibrium dynamics close to co- 
existence, we still expect that this dominant phase may 
support relatively large fluctuations of the minority (in- 
active) phase, depending on the size of the free energy 
difference and the space time "surface tension" . This was 
the scenario envisaged in [9, 10]: domains of the inactive 
phase form the slow correlated regions that form part of 
the dynamical heterogeneity in the glassy systems. 

In the soft East model, we estimate the free energy dif- 
ference between active (equilibrium) and inactive phases 
to be Afi w s*(fc(s = 0)) ~ 76 per site, per unit time. 
Both 7 = e"''/^ and e = oT^ decrease on cooling, so 
this effect promotes fluctuations of the inactive phase, 
and enhanced dynamical heterogeneity. However, we ex- 
pect the low-temperature crossover at Tx to be relevant 
for the low-temperature behavior (recall Sec. I). Assum- 
ing that (24) holds, then Cc tends to zero faster than 
e ^ c~^/'^ on cooling, with the result that e > ec at low 
temperatures, and there is no space-time phase transition 
for T < Tx 

To summarize, facilitated dynamics dominate the sys- 
tem for Tx < T < To, and this is the regime in which 
space-time phase transitions are observed. Outside this 
regime, the dynamics of the system have weaker correla- 
tions in space and time, and the response to the field s 
is smooth, not singular. The situation as a function of 
temperature T and bias s is sketched in Fig. 4. However, 
as discussed in Sec. I, the critical point in the vicinity 
of Tx is probably hard to observe in atomistic models 
of glass-formers, since supercooling the system to such a 
low temperature is likely to be very difficult. 

We end with a few comments on the space-time phase 
diagram we have proposed in Fig. 4. Firstly, while two 
critical points appear when the phase diagram is plotted 
in these variables, phase diagrams in the (e, s) plane for 
any fixed 7 are always of the form shown in Fig. 1(a). 
Except for the point where e = and 7 — >■ (the original 
East model as T — > 0), these critical points are all of 
the same type. We expect these to be in the universality 
class of the Ising model in (1 -I- 1) dimensions, as in the 
soft- FA case [20] . (Settling the universality class of these 
transitions requires an assumption that the microscopic 
left-right asymmetry of the East model is irrelevant on 
very long length scales near the critical point: we have 
not been able to demonstrate this explicitly but it seems 
plausible from our numerical results). Finally, we note 
that passage from active to inactive states and vice versa 
may be achieved by varying s and T in Fig. 4, without 
ever passing through any phase transition line. Such a 
situation requires that the symmetry properties of the 
active and inactive phases must be the same, and forbids 
any spontaneous breaking of translational invariance in 
the inactive phase. This is consistent with the behavior of 
the soft-FA model: in the analogy with atomistic glass- 
formers it requires that the inactive phase be a liquid, 
and not a crystal or glass phase where non-trivial density 
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profiles persist indefinitely. 
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